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Contact with some Machine Simuletion Models (war cames) 
Peewamaicated, that the neec for “random numbers” is real and 


pressint, yet the information aveilable to the war samers 


rh 


Gogee-ning tie Becpenties of these numbers is oO 


oO 


ten cuite 
Maguc and Limited: 

The average of the results of several runs of a war game 
is often used to evaluate the many relationships tested by 
the game. Since random numbers are the device that determines 
che yes/no answer at the many decision points of the zame, 


the faithfulness with which these numbers exhibit the prover- 


= 


ties cesired by the war gamer cirectly affect the validi 
tne Seeutits. 


In view of the critical role in war gaming played by 


these numbers it seems useful to vrovice a tested method of 


procucing the numbers, and to nrovice also the means of test- 


$e 


ng any other number-producing methods considered uscful. 
Rancom numbers have an important anplication in statis- 

tical samplines, wav caming, anc other monte carlo techniques. 

The applications mentionec in this paper will be primarily in 

tne war game area. 
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are some initiel orohleme in the study of war gaming. One 
nroblem Lies in Pavine vrecise terminology end notation to 
saplain concents that are larzelyv nrobabilistic in nature. 
Another lies in the stochastic structure of the game, i.e., 
the decision points in the game which reauire drawing a ran- 
dom number to determine the occurrence or non-occurrence of 
an event, This provices the requirement for a random number 
generator in the game. /. third oroblem is, given a random 
number generator, are the characteristics of the numbers pro- 
ducec by this generator such that the results of the game are 
valid? 

The secon: and third problem areas mentioned, that is, 
eneretion of random numbers, and testing of rancom numbers 


are the subject of the investigation in this thesis. The 


nvestigation included a search of the literature and accu- 


pds 


mulation of reference material, writing CDC 1604 FORTRAN pro- 
grams for seneration of rancom normal, exponential and chi- 
square numbers, writing programs for tests of random numbers, 
and running the progrems for the senerators and tests on the 
CPC 1604 computer. 

Sections 1 and 2 of the thesis contain cdiscussions of 


seneration of random numbers and testing of randon numbers, 


respectively. Section 3 contains the results of tests pver- 


bie 





formed on two viermai rancom number generators, an exponen- 
tial and a chi-square random number generator. Section 4 ° 
presents the conclusions reached in this paper. Appencices 
A, B and % provide justification for the methods of genera- 
tion of random numbers and Appendix P reproduces the pro- 
ecams in CDC 1604 FORTRAN Language for the four generators. 

Tt is suggested that the reader familiarize himself with 
the explanation of notation and definitions of terms immediate- 
ly following this introduction. [t is also suggested that the 


oo 


U.S. Naval Vostgraduate School thesis "Random Number Genera- 
tion on Picsital Computers" by J. ti. Barron 19 be consulted 
‘ J 


for information concerning the generation of uniform random 


numbers. 





DEFINITIONS OF TERMS AND NOTATION 
In this paper a rancom variable, X, and its associatec 
cumulative distribution function F, will be identified as 
(X,F) or, in some instances, where it is more convenient to 
refer to the associated density function f, we shall write 


(X,f£). In conformity with standard notation 


F(x) 


it 


Pr(X¥ =x) 


fx) Ps ores) 


Random Numbers: It is considered that there is available 
some formula or mechanism or device which procuces numbers 
which are regardec as drawn from a population (X%,F) or equi- 
valently (X, £). Such a number is frequently identified as a 
pseudo-random number in the Literature; however, here we shall 
identify it simply as a random number, 

We may say that a set of such random numbers is associ- 
ated with the random variable (“,F) or, more briefiy, that 
they are associated with the distribution F. 

Sequence of rancom numbers: As the formula or the generating 
device or mechanism employed produces a set of random iumbers, 
these appear serially or sequentially. We shall refer to 
these mimbers as a sequence of random ninbers, 


In referrine to random numbers, as cefined above, as- 





sociated with snecific random variables and their associated 
cumulative distribution functions, we shall need additional 
specific definitions and notation. 

NCULL,G): A normal distribution with mean [ and standard 
deviation O as parameters, 

Ule,b): A uniferm distribution over the interval [a,b]. 

In particular, U€0,1) will be used extensively. 

Randon NCL, OG) number: 4. random number associated with the 
random varieble (X,F) where F is N(lL,O). 

Random UCO,1) mumber: A rancom number associated with the 
random variable (“%,F) where F is UC0,1). 

Normal Peviate: A random N€9,1) number which has been obtained 


as a transformation of a random J(0,1) number. 


co 


senerator: «4 formula or method for procucinz a sequence of 
numbers, 

NORMGEN: The computer program to generate normal deviates 
by the Box-Muller method, 

NORMSUM: The computer program to generate normal deviates 


by the "sums of uniforms" method, 


Be eGhN: The computer program employing Marsaglia's method 


A 


to generate exponential random numbers, 
CHIGEN: The computer program to generate chi-square random 


numbers from sums of squares of independent normal deviates, 





with the number of degrees of freedom as an arbitrary 


parameter, 





Seer ron | 
GENERATION CP RANDOM NUMBERS 

Lees LC THEO 

For generation of rancom numbers associated with seneral 
rancom variables (X,F) from random UC0,1) numbers, it is nec- 
cessery to establish a relation between the uniforn distribu- 
tion and other distributions. Particularly a relation between 
the uniform and normal distributions would be useful, because 
other distributions such as the chi-square, F, and Student's 
t can be derived from the normal distribution, «A general re- 
lation between the uniform distribution and other distribu- 
tions does exist and is called the vrobability integral trans- 
formation. This relation may be stated as follows: (See Is | 
for a proof} 

Theorem: Let % be a continuous random variable 

with distribution function G, Then the randon 

variable Y ="3(X) has a density function given by 
| S Cay. 
fe 

» elsewhere; 


its distribution function is given by 


fo yo 
pee y| = ly O<y=l 
| 1 y=l 





The relation y = 3(«) is called the probability integral 
transformation, anc since %& is a non-decreasing function it 


follows that ‘see page 313 of Sie 
a - 
Pr fe oe «| = Pr LY =6(x)| = Gian, 


The above provides the basis for generating random nam- 
bers associated with an arbitrary distribution by transform- 
ing UC0O,1) numbers. 

In the case of the exponential distribution where 
mae) = 1 ave it is possible to solve explicitly for x, 
that is, x = see ley =O. In this case we say that G has 
been founc in closed form, meaning that «x has been expressed 
in terms of the logarithm function, one of the functions 
which, Like the sine and cosine, we choose to call elementary 
Pumctidonss 

For a general distribution function, it is not always 
possible to find the inverse function in closed form, that is, 
it is not always possible to exoress the inverse in terms of 
so-called elementary functions. 

For the normal distribution, a closed form expression 
for the inverse’ function is not available. In such a case 
it may be possible to employ a graphical interpretation. 


srenhical methods, however, have a minimum of decimal place 





accuracy and lead to table took-un procedures in high sveed 
digital computers, “nterins the ordinate of Figure 1 to de- 
termine abscissa values by use of a table, sometimes called 
the “Normal Curve of “rror Table", also leads to table loox- 


un methods in the computer, 


i) 


a 


Figure 4 
The disadvantage of the table look-up method of supplying 
rancom numbers for a war came is apparent when the memory 
soace Limitation of a comnuter is considered, 

ven thouch no inverse function may exist in closed 

form, it is vossible to construct generating schemes that are 
fast and have acceptable pronerties. Two such senerating 
schemes for the normal distri 


DUCIOn Wiese civen. in the next 


Geom cections. 


1,2 Normal Numbers by the Box-luller Method CNCPMG™N) 


o« 


¢ t e@ 
Bow and luller 1 proposed a methoc for generating 4 





vair of indepentent rancom N(9,1) numbers from a pair of ran- 
com U{O,1) aumbers. The transformation (see Appendix 4) for 


e ® ‘ e 
generating the pairs of (0,1? numbers is: 


ie 
Ko eee ioe 1] YS Sos 2 710 
- -e 1 2 
: T 
K = Ge@eloc U Sian Ze 
2 “e 1 2 


~4? 
dent sco 


where U Uy are indenendent JC0,1) and as ave incepen- 


In view of the four functions (sine, cosine, log. 
square root) which must be computed in the Box-iMuller method, 
this method cannot be considerec fast, although a vair of 
numbers is produced so tniat computation time must be halved 
ror comparison with other methocs. The major advantage of 
the Box-ltuller method is that it is en analvtically exact 
transformation (of pairs of nurbers) from the uniform ANeee ie 
bution and will faithfully produce a normal distribution 
eiven a iniform cistribution. 

Cne disafvantage, other than consutation time, is that 
the Sox-luller methoc will reflect whatever bias or contam- 
ination is present in the uniform cenerator. 

“hen written es a subroutine anc included in a service 


Library, thie 3o0:-tuller normal zenerator may be employed in 





two ways: The vair of NC€O,1) numbers may be delivered simul- 


ww 


eaneois ls, A in the “s-rezister and = in the %-re¢ 


the subroutine may cave “, until the next call for a randon 


ister, or 


mimber, In the first instance the user would be required to 
save K, for any case where only a single number were desired, 
In seneval, standard Library subroutines are very expen- 
Sive in machine time, as they are required to be able to ac- 
cept the ceneval cese..,.,...they must beer for sion and mag- 
nitude in many cases before beginnings the computation, In 
special cases where the ar.zument with which the subroutine is 


ao =f 
aan 


enteszed is O2%=S1, or perhaps .=277, the tests are not 
sermane and precious machine time is wasted. In cases such 


as tis (ae Euncttonal Stibroutines can be reWritten and 


}-~<4 
fe 
ty 
— 
“ 


shortened and inc part of the cenerator subroutine, 


BH - 2 @rue a “ oo = es 7: 
However there is trade-off here, as this c: 


a) 


eatly increases 


| 


the Lancet! i 


a Of fie Generator anc consequently the librery rro- 
Prem, “RA@ libracy funetion routines are then partially du5li- 


ecatec within the cenerator routine when the generator routine 


becomes wart of the library. °f the four cenerators described 


* ae = ® os oe at a aa ° = ° a 

in this thesis, these remarks nertain almost exclusively to 
a af Wien eo i peNny, 

‘heme: ane Wullor “OS UIGEN, 





1.3 Normal Numbers from Sums of Uniforms (NCR{SUM) 


i 


ev 
ne 
{fh 


cussion of the jus itaccet ton for approximation of 
a normal (0,1) random number by the sum of 12 uniform (0,1) 
random numbers appears in Appendix A, The approximation can 
be improved by increasing the number of uniform (0,1) nun- 
bers used, but machine time places a restriction on the de- 
sirability of usine more numbers. The sums of uniform random 
variables converge more rapidly to a normal distribution than 
the sums of random variables from a non-symmetric distribution. 
The proeran for cenerating these numbers is short in both 
progran steps and machine time. There are no Library sub- 
routines which must be called in, 
This method is anpealing to the intuition as there is 


the factor of "smoothing" or averaging whatever bias or con- 


tamination may be present in the uniform (0,1) generator. 


1.4 Uxnonential Random Numbers CEXPOGEN) 


marsaglia's [2 method for generating exponential ran- 
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dom numbers is more complicated in programming details than 


AIT. 


either the NOXIGTN or NGRISUIL, but Lends itself very well to 
the POXTRAN language. 
The exnected value of the n of MNarsaglia's method (sce 


Appendix 3) is approximately 1.523, so on the average only 








1,58 UC€O,1) numbers are required fron which the minimim must 
be selected, and 1.5% @iscriminations (see Appencix B) are 
sequired to assign a value ton. Assisning a value to nm also 
recuires on averace of 1.58 discriminations, Generation of 
an exnonential number thus requires an average of approxi- 
mately 4.75 UCO,1) numbers. 

This method of nroducine exnonential random numbers is 
sensitive to anv bias in the characteristics of the uniforn 


See aoe . 
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1,5 Chi-Scuare Random Numbers (CHISEN) 

Suppose that ten chi-sauare random nunbers with five 
decrees of freecon are desired and that the NCHMSUI! cenera- 
tor is use to nrocuce the normal tmumbers, It takes 12 uni- 
form mainbers tio nroduce one normal number enc five normal 
numbers to produce one chi-square number wit! Five deerecs 
of freedom. Thus we require 50 normal numbers and 600 uni- 
form numbers to nrocuce ten chi-square numbers with a rela- 
wvoemall number of dervrees af freedom, Cre has no diffi- 
culty seeing that vroduction of chi-square numbers is exnen- 


ive im machine time, ‘iowever, in war gaming the neec for 
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Vi-square numbers is Likely to he considerably less than 


mim Wee soceiini form Or nornal numbers. 





The chi-eeuare gimeretor rrocram is again relatively 
12 on ws 4 6 nit e oo pe 4 rr} 
Cirect ance no snecial cifficulties ere encountered. The 
frecuency enc coocness-cf-fit tests co present a special 
as chenging the cecrees of freedom in effect creates 
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&@& new cietribution anf’ the class intervals and expected fre- 
quencies of occurrence within the class intervals chanze. 

Zney co not chanze in such a fashion that PO-LOCTS may be 

most cconomically employer, The program becomes quite Lenethy 
in machine time if chi-scuere numbers with several cifferent 


Me q 


@@erecs Of Freeron are teste. 


1,6 Somouter Lancuage 


SOTRN was the computer Lancuage employed for all pro- 


¢ 
a 


Brans written for this thesi 


Gr 
Oe 


Pee eal i Zed at FORTRAN 


wh 
Me 


may not be the best lancuace for all tynes of war cames 


thous. successful war cames have been written in this len- 


euace, Speec anc menory space Limitations may dictate the 
use of assembly routines as onposec to FORTRAN-tyne compilers. 
Rowever, within a SI TRAN nrogvran, strines of symbolic machine 


Tae wer TY 


Mancuece may be inserted wnere the lezal FORTRAN statements 


ere not appropriate. In general, the particuler character 
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of the simulation will sugsest the most convenient lLanzuace. 


ror use in a war gane, the four random number senerators 


? q 


Poe teeeesie mist bo translatec imto the lanzuace in whieu 





the game is written. “ne word of caution is in order at this 
noint. The method used to senerate the U{0,1) numbers de- 
nends on the computer having a 48-bit word Lensth, and the 
method must be revised for use on a machine with a different 
word ieneth, 

The Control Date Corporation 1604 stored vrogram, gen- 
eval vcurpose cigital computer with a storace capacity of 


32,765 4°-bit words was available at the Naval “ostgraduate 
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School for this project. 
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nerimnents for estimating the varameters of mixec nonulations 
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4 


(3 


oh a i 9 ? a fon s 4 a y a, o o a ee ae 2 
eenerator, Tre employment consicercc in this paner will be 


an : meat s a “~ “te 
Geilo oste ws Lt tne war pane rea, 


™ 
‘ 


T£ the F with which the sequence of random numbers is 


feebe associcted ic (0,1), it is obvious that a mean near 
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zero and a veriance rear one must be demanded, but is this 
Sebi otrete Fe answer is emphatically: no, as Figure 2 
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The ensity fumetion grannec in Figure 2 could well have ¢ 


mean of sero an? @ variance of one, but the stew to the left 


root 


eliminates it from consiceration. 
test of the ceriuence of ranibers by computation of the 


first three central moments gives the criteria for the nean, 
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how we have establishec a requirement for at least com- 
nutins the first™four central moments to provide data for de- 
cidine if a sequence of numbers can be recarded as a sample 
from a cesire“ distribution, F. 

“at this point the experimenter would co well to apply 
frequency test te the sequence cf numbers by dividing the 
theoretical rance of the numbers into class intervals and 
Counting eme numbers which fall into each class interval. 
Since the theoretical frequency of occurrence of numbers in 
each class interval is known, a soodness-of-fit test Like the 
chi-square test nay be aprlied to test for a significant dif- 
Ference between samnle frequencies anc observec frequencies 
in the class intervals. 
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i ase of the N(O,1) distribution, any Sample con- 
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eidered to have been “rawn at random should have an approxi- 
mately equal number of positive numbers and negative numbers. 
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cuer.ce are generated is also of primary importance. ‘for ex- 
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The over-all effect of this sikuation is that the game must 
De played, - propjsitively Lacce number of times to obtain 
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“yo above seneral testine infornation is summarizec in 
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nronerties, even when 2éministered collectively. 
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apnearance when nlotted in 2 histogram. 
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APYENOE, oO 
EXPONENTIAL RANDOM NUMBERS 

The exponential distribution figures prominently in 
particle or radiation studies, reliability, Life testing and 
like areas which can be of interest to persons using high 
speed computers and monte carlo techniques. A power series 
expansion to compute the logarithm of a uniform random vari- 
able may be used to generate exponéneeM random variables, 
but a@ more rapid method is desired, 

G. Marsaelia of the Boeing Seientifie Research Labora- 
tories offers a simpler device for producing exponential 
rancom variables by performing discriminations on the rela- 
tive magnitudes of uniform (0,1) random variables. See D. 

The idea is to choose the minimum of a random number of 
uniform random variables, then add a random integer---say, 
let n and m be random integers according to the following 
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Let. x,,~ eo» 9%, be FY independent rancom variabies 

normally distributed, each with zero mean and unit variance. 
ue : 

Let X equal the sum of the squares of these ) rancom 

wer tallies piesa < ioe a + ay TOT a eh . X is also 

a random variable because it is a function of random vari- 


ables. The density function of the chi-square random vari- 
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has a chi-sctueve distribution with J cezrees of freedom. 
Mm example of such a random variable is as follows: 
A missile is sent toware a tarreet. If the rance and de- 
flection ecrors of the missile are independent normally 
cistributed rancon veriables each with zero mean anc equal 
variances, the squaree distance between the explosion and 
the tarzet, divided by the variance, has a chi-sauare Cis- 
tribution with two dezrees of freedom. This can be easily 
extendec to three dimensions, 
* subroutine available to a war gamer to supply inde- 
2 
nendent X numbers with arbitrary decrees of freedom thus 


can ‘ave immediate application. 
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POllowine is the Yoxw-lliller methoe of generating pairs 
cf Ni0,1) numbers from pairs of UC9,1) numbers, programmed 
in the Tontrol Pata Corporation 1604 FORTRAN lancuace. 


The program ec written merely commutes 100 numbers 
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